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Abstract Porcine deltacoronavirus (PDCoV) in Thailand 
was first detected in 2015. We performed a retrospective 
investigation of the presence of PDCoV in intestinal sam- 
ples collected from piglets with diarrhea in Thailand from 
2008 to 2015 using RT-PCR. PDCoV was found to be 
present as early as February 2013. Phylogenetic analysis 
demonstrated that all PDCoV variants from Thailand differ 
from those from other countries and belong to a novel 
group of PDCoV that is separate from the US and Chinese 
PDCoV variants. Evolutionary analysis suggested that the 
Thai PDCoV isolates probably diverged from a different 
ancestor from that of the Chinese and US PDCoV isolates 
and that this separation occurred after 1994. 


Porcine deltacoronavirus (PDCoV) causes an enteric dis- 
ease in pigs that is characterized by diarrhea and is clini- 
cally similar to porcine epidemic diarrhea (PED) and 
transmissible gastroenteritis (TGE) [4]. PDCoV, a recently 
emerged RNA virus of the family Coronaviridae, genus 
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Deltacoronavirus, was first reported in Hong Kong in 2012 
during an investigation to identify novel coronaviruses 
[14]. In February 2014, PDCoV was first detected in pigs 
with clinical diarrheal disease on farms in Ohio, USA. A 
subsequent retrospective investigation reported revealed 
that PDCoV had been present in the US as early as 2013 
[11]. PDCoV has since been identified in Canada, China, 
South Korea and Laos [1, 6, 9]. In Thailand, PDCoV was 
first reported in 2015 [10]. However, it is unknown whether 
PDCoV was present in Thailand prior to 2015. We there- 
fore conducted a retrospective study to detect the presence 
of PDCoV intestinal samples from Thailand that had been 
submitted for the identification of enteric diseases. Positive 
samples were subjected to full-length genome characteri- 
zation followed by genetic analyses to compare Thai 
strains of PDCoV to those from other countries and to 
investigate the genetic similarity and the evolutionary 
relationships among them. 

A total of 241 intestinal samples were assayed for the 
presence of PDCoV using RT-PCR. Intestinal samples 
were from 3- to 4-day-old piglets with clinical diarrheal 
disease collected from January 2008 to December 2015 and 
stored at -80 °C. The presence of other viruses causing 
enteric disease, including porcine epidemic diarrhea virus, 
and transmissible gastroenteritis virus was detected using 
RT-PCR. In brief, viral RNA was extracted from the 
samples using NucleoSpin RNAVirus (Macherey-Nagel 
Inc., PA, USA) in accordance with the manufacturer’s 
instructions. Complementary DNA (cDNA) was synthe- 
sized from the extracted RNA using M-MuLV reverse 
transcriptase (New England BioLabs Inc., Ipswich, MA, 
USA). PCR amplification was performed on cDNA using 
specific primers for the N gene of PDCoV, the spike gene 
of PEDV, and the N gene of TGEV as described previously 
[5, 10, 12]. 
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Samples positive for PDCoV were then subjected to full- 
length genome sequence characterization. Twenty-six pairs 
of primers were used to amplify the different regions of 
PDCoV [8]. PCR amplification was performed using Plat- 
inum® Tag DNA Polymerase High Fidelity (Invitrogen, CA, 
USA) in accordance with the manufacturer’s protocol. The 
PCR products were purified using a NucleoSpin Plasmid kit 
(Macherey-Nagel Inc., Bethlehem, PA, USA). The 5’- and 
3’-terminal regions were determined by using a kit for rapid 
amplification of 5’ and 3’ cDNA ends (5’ and 3’-RACE, 
Clontech, Japan). Sequencing was performed at First 
BASE Laboratory SdnBhd (Selangor, Malaysia) using an 
ABI Prism 3730XL DNA sequencer. 

The full-length PEDV genome sequence was assembled, 
and nucleotide and deduced amino acid sequences were 
aligned and analyzed using the CLUSTALW program [13]. 
The percentages of homology between the sequences at the 
nucleotide and amino acid levels were calculated. To 
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investigate its evolutionary relationships to foreign PDCoV 
strains and estimate the divergence time of PDCoV in 
Thailand, phylogenetic analysis of the full-length nucleo- 
tide sequence of the PDCoV isolate collected in this study, 
together with 21 other PDCoV isolate sequences (Supple- 
mentary Table 1), was performed using the Bayesian 
Markov chain Monte Carlo (BMCMC) method imple- 
mented in the program BEAST v1.8.3 [2, 3]. A BEAST run 
was performed based on the GIR+G+I model with a 
coalescent Bayesian skyline tree prior and a clock model 
for each analysis using 200 million generations with sam- 
pling of every 10,000 generations and the first 10% dis- 
carded as burn-in. Tracer v1.6 was used to confirm that 
post-burn-in trees yielded an effective sample size (ESS) of 
>200 for all parameters. The resulting tree was viewed and 
generated in FigTree v1.4.2. 

The results of testing for PDCoV, PEDV, and TGEV in 
intestinal samples are shown in Table |. Of 241 intestinal 


Table 1 Samples tested for porcine deltacoronavirus (PDCoV) by PCR in Thailand in 2008-2015 


PDCoV positive (%) PEDV positive (%) TGEV positive (%) 


Year Region Province Number of samples 
2015 West Rachaburi 24 
Nakhon Pathom 14 
Middle Saraburi 10 
Lopburi 6 
East Chonburi 1] 
Chachoengsao 2 
Northeast Nakhon Ratchasima 8 
Buriram 13 
2014 West Ratchaburi 2 
Nakhon Pathom 3 
Middle Saraburi 4 
East Chonburi 1 
2013 West Nakhon Pathom 4 
Ratchaburi 6 
East Chonburi 7 
2012 West Ratchaburi 4 
Middle Saraburi 2 
East Rayong 1 
2011 West Ratchaburi 68 
Nakhon Pathom 1 
Middle Saraburi 9 
East Chonburi 5 
Northeast Nakhon Ratchasima 3 
2010 West Ratchaburi 9 
Northeast Nakhon Ratchasima Z 
Udon Thani 2 
2009 1 
2008 West Nakhon Pathom 9 
Total 241 
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1 (4.17%) 22 (91.67%) 0 (0%) 
0 (0%) 10 (71.43%) 0 (0%) 
1 (10%) 10 (100%) 0 (0%) 
2 (33.33%) 6 (100%) 0 (0%) 
5 (45.45%) 11 (100%) 0 (0%) 
0 (0%) 2 (100%) 0 (0%) 
0 (0%) 3 (37.5%) 0 (0%) 
0 (0%) 13 (100%) 0 (0%) 
4 (33.33%) 6 (50%) 0 (0%) 
0 (0%) 3 (100%) 0 (0%) 
0 (0%) 4 (100%) 0 (0%) 
0 (0%) 1 (100%) 0 (0%) 
0 (0%) 4 (100%) 0 (0%) 
2 (33.33%) 4 (66.67%) 0 (0%) 
0 (0%) 7 (100%) 0 (0%) 
0 (0%) 4 (100%) 0 (0%) 
0 (0%) 2 (100%) 0 (0%) 
0 (0%) 0 (0%) 0 (0%) 
0 (0%) 9 (13.24%) 0 (0%) 
0 (0%) 0 (0%) 0 (0%) 
0 (0%) 5 (55.56%) 0 (0%) 
0 (0%) 1 (20%) 0 (0%) 
0 (0%) 1 (33.33%) 0 (0%) 
0 (0%) 6 (66.67%) 0 (0%) 
0 (0%) 1 (50%) 0 (0%) 
0 (0%) 1 (50%) 0 (0%) 
0 (0%) 1 (100%) 0 (0%) 
0 (0%) 7 (77.78%) 0 (0%) 
15 (6.22%) 144 (59.75%) 0 (0%) 


Porcine deltacoronavirus in Thailand 


Table 2 Gene lengths in four Thai PDCoV isolates and positions of amino acid deletions/insertions and substitutions in Thai PDCoV isolates 
compared to the other two PDCoV groups 


Thai isolates 


ST1_0213 


ST2_0214 


TT_1215 


NT1_1215 


Length 


China group 


US group 


Length 


China group 


US group 


Length 


China group 


US group 


Length 


China group 


US group 


Nucleotides 

Amino acids 
Deletions* 

Insertions* 

Number of substitutions 
Deletions* 

Insertions* 

Number of substitutions 
Nucleotides 

Amino acids 
Deletions* 

Insertions* 

Number of substitutions 
Deletions* 

Insertions* 

Number of substitutions 
Nucleotides 

Amino acids 
Deletions* 

Insertions* 

Number of substitutions 
Deletions* 

Insertions* 

Number of substitutions 
Nucleotides 

Amino acids 

Deletions* 

Insertions* 

Number of substitutions 
Deletions* 

Insertions* 


Number of substitutions 


* The numbers indicate the position in each gene 


© Amino acid substitutions at different positions 


Gene 

ORFla/Ib S E M N 
18,786 3,477 249 651 1,026 
6,262 1,159 83 217 342 
4017 402 ang >8pyG7© ST1yy : - : 

- SIN - - - 

91 35° Z 1 5 
4017 x 402 ang 8pyG7© S71 y7 7 - : 

96 33° : 1 5 
18,786 3,477 249 651 1,026 
6,262 1,159 83 217 342 
401) 402 ang 8pyG7© ST1yy - - ” 

7 SIN 7 - - 

91 37° : 1 5 
4017 x 402 and 8pyG7© S71 yz : - - 

96 35° : 1 5 
18,786 3,477 249 651 1,026 
6,262 1,159 83 217 342 
4017 402 ang 8pyG7© 7 ; - f 

90 D3° : 1 4 
4017 402 ang >8pyG7© - - 7 - 

95 23° Z 1 4 
18,786 3,479 249 651 1,026 
6,262 1,160 83 217 342 
4017 402 ang M8pyG7© : - - : 

- SIN - - - 

139 ae : 1 3 
4017 402 ang 8pyG7© - 7 ; - 

144 27° : 1 3 


samples tested, all were negative for TGEV. Fifteen sam- 
ples (6.22%) were positive for PDCoV. One sample that 
was collected in February 2013 was PDCoV positive. The 
sample was from a farm in Ratchaburi, a province in the 
western region of Thailand where pig density is high. 
Interestingly, the sample was also positive for PEDV. Four 
samples collected in 2014 from two farms in Ratchaburi 
were PDCoV positive. On one of these farms, samples 
from 2013 and 2014 were positive for both PDCoV and 
PEDV. In 2015, PDCoV was detected in eight samples 
collected from five farms located in Ratchaburi, Saraburi, 


Lopburi and Chonburi. These four provinces contain three 
major swine-producing areas of Thailand. These results 
suggest that PDCoV is widespread throughout Thailand. It 
is noteworthy that 144 samples (59.75%) were positive for 
PEDV and that samples that were positive for PDCoV were 
also positive for PEDV. The results suggest the presence of 
mixed infections with PEDV and PDCoV in pigs. The 
production records of the herd in which PDCoV was first 
detected showed that there had been three diarrhea out- 
breaks during 2014 and 2015, and both PEDV and PDCoV 
were detected. In three diarrhea outbreaks, the pre-weaning 
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Table 3 Pairwise nucleotide and amino acid sequence identity values for the full-length genome sequence of the Thai PDCoV isolates from 
2013, 2014 and 2015 compared to PDCoV variants from China and the USA 


Year China PDCoV variants 
Nucleotide (%) 


2013 and 2014 97.2-97.8 
2015 96.9-97.8 


mortality rate was 56.9, 52.1 and 65.3%, respectively. 
However, we were unable to define the relationship 
between PEDV and PDCoV regarding to the primary or 
secondary causes of the diarrhea outbreaks. Whether co- 
infection with these viruses would necessitate more-com- 
plex control methods or would result in more-frequent, 
repetitive outbreaks requires further investigation. 

Of the 15 PDCoV-positive samples, one from 2013, one 
from 2014, and two from 2015 were further characterized 


0.9993 


.G07 


1769 


1975 1980 1985 1990 1995 2000 2005 2010 2015 


Amino acid (%) 


93.8-95.2 
93 .2-95.0 


KP757890.1/Porcine_deltacoronavirus_isolate_CHN-AH-2004 | 
JQ065042.1/Porcine_coronavirus HKU15_ strain LHKU15-44 
KT266822.1/Porcine_deltacoronavirus_strain_CH/Sichuan/S27/2012 
KT336560.1/Porcine_deltacoronavirus_isolate CHN-HN-2014 
JQ065043.1/Porcine_coronavirus HKU15_ strain .HKU15-155 
KP757891.1/Porcine_deltacoronavirus_isolate_ CHN-HB-2014 
KP757892.1/Porcine_deltacoronavirus_ isolate CHN-JS-2014 

— KR131621. 1/Porcine_deltacoronavirus_isolate_ PDCoV/CHJXNI2/2015 


KT021234.1/Porcine_deltacoronavirus_strain_CH/SXD1/2015 “4 


KJ481931.1/Deltacoronavirus_PDCoV/USA/Illinois12 1/2014 


US PDCoV variants 


Nucleotide (%) Amino acid (%) 


97.3-97.4 
97.0-97.3 


94.2-94.4 
93 .5-94.2 


by full-length genome sequencing. The full-length genome 
sequences of the four Thai PDCoV isolates, one isolate 
detected in 2013 (P1_13_ST1_0213), one isolate detected 
in 2014 (P2_14_ST2_0214) and two isolates detected in 
2015 (P23_15_TT_1115 and P24 _15_NT1_1215), are 
25,400-25,404 nucleotides (nt) in length, with a genome 
organization similar to that of all previously reported 
PDCoV genomes: 5’UTR-ORF1la/1b-S-E-M-Nsp6-N- 
Nsp7-3’UTR [7, 14]. The ORFla/1b, S, E, M, Nsp6-N- 


P1_13_ST1_0213/PDCoV/0213/Thailand_complete_Genome 

P2_14 ST2_0214/PDCoV/0214/Thailand_complete_Genome 

P24 15 NT1_1215/PDCoV/2015/Thailand_complete_genome 
P23_15_TT_1115/PDCoV/2015/Thailand_complete genome 
KU051649/Porcine_deltacoronavirus_strain_PDCoV/Swine/Thailand/S5015L/2015_complete_genome 


P1_16 BTL_0116 


+ China-like Group 


KJ769231.1/Swine_deltacoronavirus_ OhioCVM1/2014 = 
— KJ584357.1/Porcine_coronavirus HKU15_strain_K Y4813 
KJ620016.1/Porcine_coronavirus_ HKU15_strain_M16148 


: | in 57h M82 0765.1/Porcine_deltacoronavirus_KNU14-04 


US-like Group 


KJ584356.1/Porcine_coronavirus HKU15_strain_SD3424 


KJ584358.1/Porcine_coronavirus HKU15_strain_PA3148 


KJ584359.1/Porcine_coronavirus HKU15_strain_NE3579 


Fig. 1 Phylogenetic analysis of porcine deltacoronavirus (PDCoV) based on nucleotide sequences of the full-length genome performed using the 
Bayesian Markov chain Monte Carlo (BMCMC) method. The red dots represent Thai PDCoV isolates 
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Nsp7 genes are 18,786, 3,477, 249, 651, 282, 1,026 and 600 
nucleotides in lengths, respectively (Table 2). Substitu- 
tions and deletions/insertions at the amino acid level 
occurred in several sequence regions, especially the 
ORF la/1b and S genes, of the four Thai PDCoV strains 
when compared to variants from the United States and 
China (Table 2). Interestingly, all four Thai PDCoV strains 
had amino acid deletions at *°'LK*°? and *®PVG’°". 

The pairwise nucleotide and amino acid sequence 
identity values of the full-length genome sequences of the 
four Thai PDCoV isolates are shown in Table 3. These 
sequences were 98.8-99.1% identical at the nucleotide and 
97.4-97.9% identical at the amino acid level and were more 
closely related to each other and to variants from China and 
the United States than to Thai strains isolated in 2015, 
suggesting that the Thai PDCoV population is continuously 
changing. 

Phylogenetic analysis based on the full-length genome 
sequence demonstrated that PDCoV isolates are divided into 
two groups, the Chinese- and US-like-groups, as was 
reported previously [1]. All four Thai PDCoV isolates were 
grouped in a novel cluster separate from variants from the 
United States and China but were more closely related to the 
Chinese isolates than to the US PDCoV isolates (Fig. 1). The 
currently circulating isolates of PDCoV share a common 
ancestor that was present around 1994, with 95% HPD values 
covering 30 years from 1975 to 2004 (Fig. 1). Since then, 
Thai and US/China PDCoV clades diverged from the com- 
mon ancestor and evolved independently. The emergence of 
the Thai lineage was estimated to have occurred between 
2004 and 2012, while the US/China lineage evolved earlier, 
between 1994 and 2003. The clustering of the PDCoV iso- 
lates in the Thai PDCoV clade is well supported, with a 
posterior probability of 0.9993. Thai PDCoV sequences 
isolated in 2013 and 2014 grouped separately from Thai 
PDCoV sequences isolated in 2015. Moreover, all four Thai 
PDCoV isolates were descended from a common ancestor 
between 2004 and 2014 (Fig. 1). The analysis of amino acid 
changes also suggests that they are from a different lineage. 

In conclusion, we retrospectively investigated the presence 
of PDCoV in Thailand and found that PDCoV has been present 
in Thailand since at least February 2013 and that it is present in 
mixed infections with PEDV. Interestingly, the PDCoV vari- 
ants in Thailand are different from variants in other countries 
and belong to a novel group of PDCoV strains. Evolutionary 
analysis suggested that the Thai PDCoV isolates likely 
diverged from a different ancestor from that of the Chinese and 
US PDCoV isolates, and the separation occurred after 1994. 
This study provides interesting results associated with the 
emergence and the evolution of PDCoV. Further investigation, 
including more retrospective studies of available samples, will 
provide a better understanding of this pathogen. 
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